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[57] ABSTRACT 

Techniques of combining separate but correlated measure- 
ments to form a second-order or higher order correlation 
function to suppress the effects of noise in the initial 
condition of a system capable of retaining memory of an 
initial state of the system with a characteristic relaxation 
time. At least two separate measurements are obtained from 
the system. The temporal separation between the two sepa- 
rate measurements is preferably comparable to or less than 
the characteristic relaxation time and is adjusted to allow for 
a correlation between two measurements. 

9 Claims, 9 Drawing Sheets 
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METHOD FOR SUPPRESSING NOISE IN 
MEASUREMENTS 

RELATED PATENT APPLICATION 

This application claims the benefit of the U.S. Provisional 
Application No. 60/050,701, filed on Jun. 24, 1997, the 
disclosure of which is incorporated herein by reference. 

U.S. Government may have certain rights in this invention 
pursuant to Nasa contract number NAS7-1407. 

FIELD OF THE INVENTION 

This invention relates to measurements in general and 
more specifically to methods of using correlated measure- 
ments for improving sensitivity, resolution, and information 
content of measurements. 

BACKGROUND 

The accuracy and precision in measurement are usually 
limited by a variety of noise sources. Noise may be caused 
by sources that are internal or external to a system of 
interest. For example, instrument noise and environmental 
noise are external noise. External noise can be reduced by, 
for example, improving the instrument and isolating the 
influence of the environment. Internal noise is caused by 
inherent processes within the system under test. Thermal 
fluctuations, caused by non-zero operating temperature, and 
quantum fluctuations, for example, contribute to the internal 
noise and introduce uncertainty in measured quantities. Such 
internal noise may be small compared to the external noise 
in many applications. As the instrument technology and 
measuring techniques improve, the external noise can some- 
times be reduced to a level comparable to or even less than 
the internal noise. Therefore, in certain applications, the 
internal noise may become a primary source of noise. 

The interaction of an instrument and a sample under 
measurement can also introduce noise to a measurement. 
This is called a back action effect. The accuracy or precision 
of signal detection can be degraded by this back action 
effect. The measurement problems associated with the inter- 
nal noise and back action effect can become more significant 
when the dimension of a system under test or its component 
parts are microscopic. An example is the damping of the 
motions of trapped ions in sensitive mass spectrometry 
measurements by the detection circuit used to detect those 
motions. See, Lowell S. Brown and Gerald Gabrielse, Geo- 
nium Theory: Physics of a Single Electron or Ion in a 
Penning Trap, Rev. Mod. Phys. 58 (1), 233-311 (1986). This 
damping may reduce the resolution of the mass spectrum 
and degrade the precision of the measurement. 

The effects of these noise sources in general are present in 
any system of an arbitrary scale that is sufficiently isolated 
from its environment so that other noise sources are less 
important. For devices of a molecular size, the internal nose 
and the back action effect may be more significant than 
observation and operation in devices at larger size scales. 
Noise can be suppressed in a measurement by correlation 
techniques. A first-order correlation function, for example, is 
an average of a detected signal such as an electrical current. 
Since noise is usually random, the first-order correlation 
function can be used to reduce the effect of the noise. If the 
signal has a non-zero average, the signal- to -noise ratio can 
be improved by averaging repetitive measurements of a 
first-order correlation. 

A second-order or higher-order correlation function is an 
average of a multiplication of two or more measurements of 
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a detected signal. See, e.g., Weissbluth, Photon-Atom 
Interactions , Academic Press, pp. 276-286, San Diego, 
Calif. (1989). Such a correlation function may be desirable 
when characteristic fluctuations in the signal as a function 
5 of, for example, time or spatial position provide useful 
information about the system. Given two measurements of 
a signal S b and S a made at different “positions” q in phase 
space and/or different times, the second-order correlation 
function is 

10 <S2>=<Sb(<lb,tb)Sa(<Ia,ta)>, 

where the bracket indicates an average over repetitions of 
the measurements. Such a second-order correlation function 
is of particular importance when the average of a signal over 
15 repetitions of the individual factors (i.e., the first-order 
correlation function) would not depend on either position or 
time. Higher-order correlation functions may also be used 
for measurements. 

When the two measurements differ only in the times at 
20 which they were made, the function is an autocorrelation 
function. If they are measurements of the same signed 
quantity measured at two times, then the Fourier transform 
with respect to the time difference (t fl -t fc ) is the power 
spectrum (or spectral density function) of this quantity. The 
25 connection of these quantities as a Fourier pair is the 
Wiener-Khinchin theorem. If <S 2 > is independent of the 
time of the initial measurement, the process is said to be 
stationary, such as a measurement made with the observed 
system in a steady-state interaction with the measurement 
30 apparatus and its environment. When the position and/or the 
measured quantity differ between the two measurements, the 
function is called a cross-correlation function. 

SUMMARY 

35 The present disclosure describes techniques of combining 
separate but correlated measurements to form a second- 
order or higher order correlation function to suppress the 
effects of noise in the initial condition of a target system 
capable of retaining memory of an initial state with a 
40 characteristic relaxation time. At least two separate mea- 
surements are obtained from the system. The temporal 
separation between the two separate measurements is pref- 
erably comparable to or less than the characteristic relax- 
ation time and is adjusted to allow for a correlation between 
45 two measurements. 

One embodiment decomposes a measurement into at least 
four distinct epochs, where all of the epochs, except possibly 
the first, are comparable to or less than any relaxation time 
of the system under measurement. In a first epoch, typically 
50 no external probe or excitation is applied to the system. In 
a second epoch, a first measurement of some observable of 
the system is performed. In a third epoch, the system evolves 
with or without an external perturbation. In a fourth epoch, 
a second measurement is performed to measure a second 
55 response. Next, a type of second-order correlation function 
is formed by using the first and second responses and is 
configured to suppress noise in the initial condition of the 
system. A third-order or higher-order correlation function 
may be similarly constructed. 

60 A scalar transformation may be applied to one or more of 
the responses measured in the second and fourth epochs to 
enhance suppression of noise. One implementation trans- 
forms the first response by a “dispersive” function and forms 
the second-order correlation by multiplying the transformed 
65 first response directly with the second response. 

These and other aspects and advantages of the present 
invention will become more apparent in light of the follow- 
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ing detailed description, the accompanying drawings, and 
the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1A is a diagram illustrating four epochs within one 5 
cycle of a measurement procedure according to one embodi- 
ment of the invention. 

FIG. IB is a flowchart of forming a second-order corre- 
lation function based on transformations of two different 10 
measurements made in one measurement cycle in FIG. 1A. 

FIG. 2 is a block diagram illustrating the multiplicity of 
possible sample and detector arrangements, as well as the 
multiplicity of possible detector coupling arrangements. 

FIG. 3 is a diagram illustrating an exemplary timing 15 
diagram of the measuring technique of FIG. 1A for a 
force -detected NMR measurement. 

FIGS. 4A and 4B are diagrams depicting calculated 
signals from the two measurements as functions of time for, 
respectively, two different total numbers N=10 1:L and N=10 2 20 
of isochronous spins with polarization P=10 -5 . 

FIGS. 5 A and 5B respectively show calculated detection 
results of using a conventional first-order correlation tech- 
nique and a dispersive -operator second-order correlation ^ 
technique for different total numbers of spins N. 

FIG. 6 is a diagram showing Monte Carlo simulations of 
the signal-to -quantum-noise ratio of the integrated NMR 
spectrum, neglecting relaxation, for various detection 
schemes. 30 

FIG. 7 is a diagram showing a series of Monte Carlo 
simulations of SNR vs. number of spins N using the 
dispersive -operator version of the method with uncorrelated 
instrument noise added to each measurement in an NMR 
experiment for a spin polarization of P=10 -5 . 35 

FIGS. 8 A and 8 B show approximated quantum-noise - 
limited spectra for (a) first-order pointwise NMR and (b) 
second-order pointwise NMR, respectively. 

FIGS. 9Aand 9B are flow charts for one operation process 
in implementing the present technique using the BOOMER- 40 
ANG scheme for force-detection of NMR. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

45 

The measurement techniques disclosed herein are appli- 
cable to a system on a time scale over which at least partial 
memory is retained of an initial state of the system. The 
system’s loss of this information is characterized by a 
relaxation time. At least some effect of a system’s initial 50 
state remains in the system in a time comparable to or less 
than the characteristic relaxation time, even when one or 
more perturbations are made to the system during that time. 

An ensemble of nuclear spins in nuclear magnetic reso- 
nance measurements is an example of such a system with 55 
memory. The spins have a longitudinal relaxation time T ± to 
describe the population decay from an excited state (one 
Zeeman state) to an equilibrium state (another Zeeman state) 
at a rate of T 1 ~ 1 . Hence, the deviation of a population of this 
state from its equilibrium value is substantially depleted 60 
after a time much longer than T ± so that the system no longer 
has a memory of that particular value of the population of 
the excited state. This spin system may also have a trans- 
verse relaxation time, T 2 , to describe the coherence time of 
the relative phases of the nuclear spin states. The charac- 65 
teristic relaxation time for the transverse magnetization of 
such a spin system is therefore indicated by the shorter time 
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of the T-l and T 2 . For other observables there are similar 
relaxation times. 

Other atomic and molecular systems can also represent a 
system with memory. The relevant energy states may be 
Zeeman states, atomic states, or molecular states. They 
could as well be states of oscillators, such as mechanical 
oscillators or electromagnetic fields described as quantum 
mechanical systems or any other states of a system of 
interest. 

FIG. 1A shows a time line for one repetition of the events 
according to one embodiment. A complete experiment may 
require repeating the steps indicated in the time line with or 
without modifications to these steps. When no modifications 
are made, this will be referred to as a repetition. The time 
line is divided into four epochs: a waiting period, a period 
for performing a first measurement on the ensemble, a 
general evolution period, and a period for performing a 
second measurement. 

A first epoch is a waiting period in which the ensemble is 
prepared for measurements. Examples of this epoch include 
waiting for the equilibrium nuclear magnetization of a 
sample in a magnetic field to develop or preparing a non- 
equilibrium spin state in some way. The waiting period may 
involve active external intervention to speed up the rate at 
which the state of the ensemble becomes uncorrelated with 
the state in which the ensemble is left by a previous 
repetition. There is no strict requirement that the state of the 
ensemble be fully uncorrelated with the previous repetition, 
but inclusion of such correlation in the analysis would 
generally be necessary unless it is made small by the waiting 
period. In the quantitative examples, the waiting period will 
be assumed to have eliminated correlation between repeti- 
tions. It is not required to have some initial order and/or 
some coherent excitation, so that this waiting period can be 
relatively short and can leave the ensemble arbitrarily dis- 
ordered and/or incoherent. For example, NMR measure- 
ments can be performed on systems with a high spin 
temperature at which spins become unpolarized. This is an 
advantage because creating a reproducible spin polarization, 
(for example by allowing the spin system to return to the 
sample temperature) can be the most time-consuming part of 
an NMR measurement. Especially at cryogenic 
temperatures, measurements based on spin-lattice relaxation 
may require too much time to be practical. Creating a 
disordered state suitable for the next repetition can be a 
relatively fast process. 

The second and fourth epochs are detection or measure- 
ment periods. It is during these periods that the system 
interacts substantially with the measurement apparatus and 
the response of the interaction is recorded. The system may 
be electromagnetically manipulated in a desired way, e.g., 
with incident radiation, to facilitate the measurement. Each 
measurement period can include measurements at one or 
more times of one or more quantities of interest. The second 
measurement period can be either similar to or different from 
the first measurement period. 

The third epoch between the first and second measure- 
ment periods is a period of generalized evolution. The 
generalized evolution period may contain multiple distinct 
steps. The state of the system changes either passively under 
internal interactions or actively under a desired external 
manipulation in this generalized evolution period. Certain 
parameters of the system, such as frequencies and decay 
rates, are “encoded” in the system by this evolution. 

Upon obtaining the first and second measurements, a 
second-order correlation function S 2 can be constructed. In 
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general, each of the two measurements may undergo a 
transformation prior to their being combined to enhance 
noise suppression by the second-order correlation function. 
Assume a transformation for the first measurement is rep- 
resented by a scalar function and a transformation for the 5 
second measurement is represented by f a , then the trans- 
formed measurements are f b ( S fc ) and .f fl (S a ), respectively, 
where S b and S a , respectively, represent first and second 
measurements (or the “before” and “after” measurements). 
The second-order correlation function is given by <S 2 >=<f b 
(S*)m)>. This process is illustrated by a flowchart in FIG. 

IB. 

Well known methods such as filtering and fitting proce- 
dures can be used for each of the two separate measurements 
in order to obtain the best statistical estimate of S b and S fl , 15 
given the data from the two measurement periods. The 
duration of these measurements is typically shorter than the 
measurement periods of many conventional measuring tech- 
niques. This is because of the need to complete the time line 
in a time short or comparable to the relaxation time and/or 20 
because longer measurement periods would not improve the 
sensitivity of the second order correlation function. Numeri- 
cal examples follow. Since the two separate measurements 
are made in a period during which the system still retains 
memory of the initial state, the two measurements are 25 
correlated by the initial state. This converts the random 
initial condition noise in many conventional measuring 
techniques into correlated signals that do not adversely 
affect the detection or do so to a lesser extent than in 
conventional methods. 30 

The first and second transformations f b and f a may be in 
various forms depending on the properties of the system 
under measurement. The exact transformations should be 
configured to enhance or optimize the noise suppression of 
the resultant second-order correlation function. The simplest 35 
transformations f b and f a are unity transformations, i.e., 
and ,f a (S a )=S a so that <S 2 >=<S fc S fl >. 

The four epochs shown in FIG. 1A complete one mea- 
surement cycle for obtaining two separate measurements of 
a system. One measurement cycle can be used to form a 40 
second-order correlation function. In order to achieve cer- 
tain correlation between outcomes of first and second 
measurements, it is desirable to set the temporal separation 
of the two measurements (i.e., the duration of the third 
epoch) within a range that is comparable to or less than the 45 
characteristic relaxation time so that the system still retains 
memory of the first measurement when the second measure- 
ment is performed. In addition, depending on the properties 
of the system and external manipulation in the third epoch, 
the temporal separation can be adjusted within the above 50 
range, i.e., within the range for which the two measurements 
are correlated with each other. This correlation may be 
modulated as a function of non-identical repetitions of the 
experiment and even cross through zero (i.e., become tran- 
siently uncorrelated for certain generalized evolution peri- 55 
ods that appear in a set of related measurements). In 
addition, additional pairs of the third and fourth types of 
epoch could be appended to obtain a third order or a higher 
order correlation function. 

The technique shown in FIGS. 1A and IB in general 60 
covers any and all repeated measurements of a system, with 
a period of coherent system evolution separating at least two 
of the measurements. During this period, the random devia- 
tion of the signal for each measurement from its expected 
value is at least partially correlated among the 65 
measurements, with the possibility that this correlation is 
modulated during the evolution period. The measurements 
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can be performed in different ways. For example, a single 
signal detector can be used to obtain the measurements of a 
single sample at different times; two or more different or 
similar detectors may be disposed relative to the single 
sample at different locations to obtain measurements at the 
same time or different times; or measurements can be done 
on two or more samples with one or more detectors. 

FIG. 2 is a generalized schematic of the multiplicity of 
possible sample and detector arrangements. It may be advan- 
tageous to use certain mathematical functions of each mea- 
surement (with the functions not necessarily the same for 
different measurements) as factors in a product involving 
either all of the measurements or subsets of the measure- 
ments. This product is treated as the signal that is averaged 
over for repetitions of the experiment. The correlation of 
initial-condition noise among the measurements can 
improve the signal-to -noise ratio compared to first-order 
techniques. Thus it is useful to take the products of mea- 
surements for which the underlying observable of interest is 
correlated in the functions so multiplied. FIG. 2 shows one 
embodiment having one or more independent samples and 
one or more detectors for each sample for parallel and/or 
sequential measurements. Some of the advantages accrue 
because the use of correlated measurements allows each 
sample to be made small enough that the simultaneous 
observation of many samples is practical and recommended 
relative to the alternative of combining the samples in a 
larger apparatus. 

Therefore, the technique of FIG. 1A may use multiple 
detectors which are either always coupled, conditionally 
(dynamically) coupled, always uncoupled, or combinations 
of these arrangements. The scheme may also use different 
detectors to make measurements on different parts of the 
same sample or on different samples. The scheme could also 
generate and exploit multiple aspects of each measurement, 
such as the various Fourier components of each measure- 
ment. 

One aspect of the technique shown in FIG. 1A is to 
explore the relation of certain sources of noise during a 
period in which a system retains memory. Some of these 
sources of noise have been largely neglected in many 
conventional measuring techniques. One of these is initial 
condition noise, which is defined as uncertainty in a desired 
observable which is a result of intrinsic uncertainty in a set 
of nominally perfect measurements on the ensemble of 
systems under test. Such uncertainty may include thermal 
fluctuations, which are an aspect of any system at a nonzero 
temperature, and quantum mechanical uncertainty in mea- 
surements of the system, which can persist even in the 
hypothetical limit of zero temperature. In addition there may 
be additional initial condition uncertainty due to inadequate 
control of any thermodynamic constraint, including tem- 
perature. Number-density fluctuations of molecules under 
study also contribute the uncertainty. The fluctuations in an 
initial condition that persists at fixed thermodynamic con- 
straints are equilibrium properties (the area of the spectral 
density function or the zero time point of the corresponding 
correlation function). When these fluctuations for an observ- 
able of interest are larger than the mean of that observable, 
the SNR of the first-order measurement can fall below unity 
and the fluctuations dominate the second order correlation 
function of that observable. One feature of the present 
method is a prescription for making novel time-domain 
measurements based on extensions of the second-order 
correlation function which provide better and/or different 
information, especially in this limit where fluctuations domi- 
nate mean values. 
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In order to make effective measurements, it is often 
desirable to modulate the observable of interest in some 
coherent fashion during the time in which it drives a 
detection apparatus. This has the effect, for the duration of 
the modulation, of convoluting the spectral density of the 
detector with the spectrum of the modulation. Modulation is 
particularly important for observables whose spectral den- 
sity would otherwise peak at the zero frequency, since the 
effects of numerous noise mechanisms increase at low 
frequencies. Thus, quasi-constants of the system are usually 
measured at the modulation frequency. Nevertheless, it is 
possible to use equilibrium statistical mechanics to estimate 
their distributions. Some of the preferred embodiments 
make use of modulation schemes of this sort. 

In other cases, the fluctuation of interest will, even in the 
absence of any experimental modulation of the observable, 
be at a non-zero frequency. Then the detection apparatus will 
measure that frequency range of the spectral density that 
couples to it. The second-order, SQUID NQR experiment by 
Sleator [Sleator, 1987 #43] is such a case. The steady-state 
spectral density of the rf magnetization of the system, 
unperturbed by incident rf, provided the signal. 

The measuring technique shown in FIGS. 1A and IB can 
be illustrated by measurements in force -detected nuclear 
magnetic resonance (“NMR”). Nuclear magnetic resonance 
is a known technique for detecting certain properties of 
materials having atoms with nuclear magnetic moments 
associated with non-zero nuclear spins. Nuclear spins are 
intrinsic properties of atomic nuclei and can be used as part 
of noninvasive probes for analyzing many materials. NMR 
spectroscopy and NMR imaging have been widely used to 
analyze the electronic and molecular structure, motion and 
material properties of a wide variety of samples. 

In the absence of a magnetic field, the equilibrium orien- 
tations of nuclear moments are random and the energies 
associated with different orientations of a nuclear moment 
are small and largely indistinguishable from one another. In 
the presence of an external static magnetic field, these 
nuclear magnetic moments acquire a net (average) orienta- 
tion and assume certain allowed quantized energy levels 
known as nuclear Zeeman levels. A transition may be 
induced between an appropriate pair of Zeeman levels by 
illuminating a sample with electromagnetic radiation in 
resonance with the transition between two Zeeman levels. 
The populations and/or phase coherence between nuclear 
spin states can be changed by such a resonant excitation. In 
practical applications, the magnitude of the static magnetic 
field is usually within a range such that the energy separation 
between two adjacent nuclear Zeeman levels corresponds to 
the radio frequency (“rf”) spectrum (~10 6 -10 9 Hz). 
Accordingly, a time-varying magnetic field for modulating 
the magnetization of a sample is within the rf range. 

Leskowitz et al. developed a “BOOMERANG” method of 
force -detected NMR in which a sample is placed in a 
homogeneous static magnetic field and an NMR signal is 
detected by measuring the motion of one or more magnetic 
sensors. The BOOMERANG method provides for better 
observation of magnetization, enhanced resolution and no 
gradient. See, U.S. patent application Ser. No. 08/872,528, 
filed on Jun. 11, 1997 and Leskowitz et al. “Force-detected 
magnetic resonance without field gradients,” Solid-State 
Nuclear Magnetic Resonance, Vol. 11(1,2), pp. 73-86 
(1998), which are incorporated herein by reference. In the 
BOOMERANG method, each measurement separately is 
proportional to the longitudinal spin angular momentum l z 
of some particular magnetic isotope in the sample, which is 
in a static magnetic field along the z-axis. Cyclic inversion 
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of the longitudinal spin magnetization drives a change in the 
amplitude of a nearby magnetic oscillator, and this change is 
measured. 

Two separate measurements are needed in order to obtain 
5 a second-order correlation function of the NMR signal. The 
timing diagram based on the measuring technique of FIG. 
1A is shown in FIG. 3. A first cyclic inversion is performed 
for a duration t b after a waiting period. After an evolution 
period t 1? a second cyclic inversion is performed for a 
10 duration t a . Two external spin manipulations are applied to 
the sample during the evolution period in the simplest 
example of the BOOMERANG experiment. First, a jt/2 rf 
pulse is applied to the sample to convert initial longitudinal 
magnetization into transverse magnetization. After time L at 
15 the end of the evolution period, another jt/2 rf pulse is 
applied to the sample to convert a transverse component of 
the magnetization that survives the evolution period back 
into longitudinal magnetization. The evolution period L is 
less than or comparable to the spin relaxation times which 
20 damp the coherent spin evolution induced during t 1 . One 
difference between the time line of this second-order BOO- 
MERANG and the first-order BOOMERANG experiment is 
that a driving/detection step is performed first, prior to the 
generalized evolution period t 1 . 

25 The two signals that are measured before and after the 
evolution period are fit to give, respectively, random devi- 
ates [IJ fc and of the spin angular momentum pro- 

jection at the end of the first detection period and the 
beginning of the second detection period. These quantities 
30 exemplify the quantities S b and S a , respectively. We are 
introducing here a general notation of the form [X(t)],-, 
where X is a quantum-mechanical operator and the subscript 
i includes any labels other than times per se needed to label 
the measurement. The time variables denoted in parentheses, 
35 as arguments of the operator, will be viewed as “dates”, in 
contrast to continuous variables. Their precise meaning will 
depend on the representation (e.g. Schrodinger vs. Heisen- 
berg representation) being used for calculation. The quan- 
tities [X(t)],. so denoted are, in general, obtained as best fit 
40 parameters to data obtained by a measurement that takes 
place over a duration, which for compactness is not explic- 
itly included in the notation. These durations are, for 
example, the lengths t b and t a of the respective measurement 
periods in FIGS. 1A and 3. The absence of a time label or 
45 the time zero indicates the end point of the before measure- 
ment. 

Note that the random deviates of the spin angular momen- 
tum are signed quantities. The product of these two 
estimates, obtained for each repetition of the time line, is the 
50 random deviate [S 2 (t 1 )]=[I^] fc [I^(t 1 )] fl of the desired correla- 
tion function. The expectation value <S 2 (t 1 )> of this product 
can be non-zero, even when the expectation value over 
repetitions of the separate signals is zero, as for an unpo- 
larized spin system. Thus it may be averaged over repeti- 
55 tions of the procedure to improve the signal-to-noise ratio. 

A special case is where all products of measurements, for 
a certain value of t l9 may be added together to provide an 
experimental estimate of <S 2 (t 1 )>. This average may be 
useful even when the sensitivity of the individual measure- 
60 ments is inadequate to usefully distinguish the individual 
eigenstates of the observable or to make any useful binning 
(vide infra) of the individual measurements. Repetition of 
the experiment leads eventually to usable sensitivity for 
< S 2 (t 1 ) > . The theoretical interpretation for the case where, 
65 regardless of the measured value, the products are just added 
together will be notated for a system whose state subsequent 
to the waiting period is described by the density operator 
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p(0). If the commutator [p(0), IJ=0, then <S 2 (t 1 )>=Tr{p(0) 
IXOi)} is the quantity to be calculated from theory for 
comparison to experiment. Here the operator I^(t 1 ) is the 
operator \ z propagated through the evolution period t l9 in the 
sense of the Heisenberg equation of motion, including 5 
possible dissipative effects. This formalism is used here to 
quantify the SNR advantages of the present method. 

Certain advantages and novelty of the present measure- 
ment technique can be demonstrated by using calculations of 
statistical moments of expected nuclear magnetic resonance 10 
signals for a system of N spins- 1 /^ in the BOOMERANG 
force -detected NMR experiment. This method allows for 
high-resolution spectra to be encoded during evolution. For 
this numerical example the N spins are assumed to be 
isochronous. After the possible waiting period, the first (or 15 
“before”) measurement samples a random deviate [l z \ b from 
the underlying distribution of possible initial conditions for 
the rest of the experiment. A resonant rf pulse is then applied 
to the spin system to generate coherence, followed by a 
period of coherent evolution (assumed for the present cal- 20 
culation to be free of dissipation). At the end of the evolution 
period, another resonant rf pulse is delivered to convert the 
encoded coherence into encoded spin population differ- 
ences. This step is followed by the second measurement of 
longitudinal magnetization [IJ fl . Even when the spin noise 25 
dominates the expected signal for either or both of the 
individual measurements, the spin noise will be correlated 
between the first and second measurements. This correlation 
will be modulated during the evolution period and so 
encoded with the frequencies of interest. 30 

More specifically, the Hamiltonian for the system can be 
written as 

H=-h(%I z , 

35 

where co z is the Larmor precession frequency of the spins. 
The expected value for the first measurement at thermal 
equilibrium is 

PN 40 


10 

The calculation indicates the degree of correlation 
between each “before” measurement and its corresponding 
“after” measurement. The spin noise is highly correlated in 
the vicinity of local maxima in the signal where the evolu- 
tion operator has cyclically returned to the identity operator 
(neglecting relaxation). Likewise, the spin noise is highly 
anti-correlated in the vicinity of local minima. At nodes in 
the signal, the spin noise is uncorrelated since a new 
distribution of total spin eigenvalues for P=0 is sampled at 
these t ± points by the second measurement, regardless of the 
longitudinal polarization at the start. 

The simplest procedure for exploiting the correlation in 
spin noise between the first and second measurements is to 
multiply the two measurements for each t 1 point. Referring 
to FIG. IB, this corresponds to unity transformations of both 
measurements: jfc(S fc )=S fc and .f fl (S a )=S a . Accordingly, a 
respective second-order correlation function, <S 2 (t 1 )> can 
be expressed as 

1 R 

(Site)) = -Y J Vzm„,Mh)] ll s 

K r=l 

The subscript r indexes repetitions of the experiment for a 
large total number of repetitions R. The quantum calculation 
of the corresponding expectation value of this product of 
first and second measurements is written in the Heisenberg 
representation as 

<S 2 (t 1 )>=Tr{p(0)I z (0)I z (t 1 )}=<I z (0)I z (t 1 )>. 

The quantum uncertainty in the product of the two measure- 
ments is calculated as 


AS 2 (r,) = V<[^(0)/z(fi)] 2 >-(/z(0)/z(fi)> 2 • 

The expectation value and quantum uncertainty for this 
example of the second-order method can be derived exactly 
as follows for arbitrary P and N: 


(S 2 (h)) = </ z (0)/ z (ri)> = ^ [1 - P*(N - 1)] cos^ 

AS 2 (h) = - Jn 2 + P 2 N 2 (N - 1) 4 - [N(N - 2) + P 2 (3N 3 -11 N 2 + 8AQ- P*(4N 3 -10 N 2 + 6N )] cos 2 o> L t x 
4 


where P is the spin polarization defined by P=tanh 
(hoo z /2kT). The quantum uncertainty in this measurement is 

1 f T 55 

A/ Z = -VW-P 2 ). 

FIGS. 4A and 4B depict calculated signals from the two 
measurements as a function of time, respectively, for two 
different total numbers N=10 1:L and N=10 2 of isochronous 60 
spins with polarization P=10 -5 . Open circles represent 
“before” measurements, and closed circles represent “after” 
measurements for each t ± point. The amplitudes are scaled 
independently for the two values of N. Instrument noise is 
taken to be negligible compared to the quantum noise in 65 
these calculations. The effect of non-negligible instrument 
noise will be discussed below. 


These two quantities are essentially independent of polar- 
ization in the regime where P^N«1. In this regime the 
signal-to -quantum-noise ratio for the second order scheme is 
approximated by 

V N cosw/ 1\ 

SNR q} , 2 - 

^ N + (N - 2) cos 2 w//i 

This SNR is also largely independent of N, since N appears 
to the same order in both the numerator and the denominator. 
One feature of this second-order correlation technique is that 
the quantum noise retains a L dependence even for vanish- 
ing P and large N. For the corresponding first-order “ordi- 
nary” NMR experiment, the spin noise in this regime leads 
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to SNR«1. Hence, the improvement in SNR by using the 
procedure of the present invention is on the order of (P 

VN)- 1 »1. 

Spin noise begins to dominate the ordinary signal when 
the distribution of total spin angular momentum in the initial 5 
condition starts to significantly overlap with zero. This leads 
to an increasing number of measurements in the signal 
average (ergodic ensemble average) whose sign is opposite 
to the sign of the average. The multiplication scheme 
corrects for this sign fluctuation and converts repetitions of 10 
the experiment from destructive to constructive contributors 
to the signal average. 

Other non-unity transformations and f a may also be 
used to construct various second-order correlation functions 
based on the same measurements as in FIG. 3. For example, 15 
a possible correction for the fluctuations is to multiply each 
“after” measurement for a given B point by the reciprocal of 
its corresponding “before” measurement, i.e., ^(S & )=l/S & 
and ^(S a )=S a . This division scheme has value for t ± points 
where the expected value of the “before” measurement is 20 
large compared to the total fluctuations. The quotient, 
however, suffers large fluctuations at t ± points where the spin 
noise is uncorrelated and the total fluctuations are larger than 
the average signal. 

A scheme can be used to correlation functions that com- 25 
bine the advantages of the multiplication and division 
schemes while suppressing their disadvantages. It has been 
shown using calculations that the “before” measurement 
may be transformed by a “dispersive” function to signifi- 
cantly improve the signal detection in some applications, 30 
such as BOOMERANG force -detected NMR. One embodi- 
ment of the operator for transforming the before measure- 
ment in the “dispersive-operator scheme” is 


MU) = 


Z h 

V — 

h z 


Zh 

Z 2 + 1\ 


35 


where £ is an optimization scaling parameter, which can be 40 
adjusted to optimize the SNR of the correlation function. In 
the regime where quantum noise dominates, the optimum 
choice for t, is roughly equal to the quantum uncertainty in 

I,: 

45 

<r h = 


The Heisenberg representation operator for the “after” mea- 
surement of I z is 

I z (ti)=I z cos (DjA+I* sin (o^A 

The final signal is computed as 


<S 2 (fi» = ( 


?/z( 0) 

<r 2 + PA 0) 



55 


FIGS. 5 A and 5B respectively show calculated detection 
results of using (a) the first-order and (b) the dispersive- 60 
operator scheme for different total numbers of spins N. The 
polarization is fixed at P=10 -5 for all simulations and each 
transient is scaled independently. Modulated signal is 
noticeable above the noise in the dispersive-operator scheme 
(FIG. 5B) even for a low spin density of N=100, while the 65 
first-order correlation technique of an NMR signal in FIG. 

5A requires roughly N=10 10 spins to achieve the same 


signal-to-noise ratio. This indicates an improvement of 
about eight orders of magnitude in the sensitivity of the 
present method to the prior art method under these condi- 
tions. 

FIG. 6 shows Monte Carlo simulations of the signal-to- 
quantum -noise ratio of the integrated NMR spectrum, 
neglecting relaxation, for detection schemes using the con- 
ventional first-order correlation technique, the direct multi- 
plication second-order-correlation version of the present 
invention, the dispersive-operator version, and the division 
version. The simulations assume that instrument noise is 
negligible compared to quantum noise. The improvement in 
SNR is by orders of magnitude for the multiplication and 
dispersive -operator schemes for PVN«1. Therefore, SNR 
can be significantly improved over the first-order correlation 
technique by using properly constructed second-order or 
higher order correlation functions. 

One aspect of the invention is to use properly-formed 
second or higher order correlation functions to improve the 
SNR of measurements even when instrument noise domi- 
nates the internal noise (such as quantum noise). The instru- 
ment noise need not be correlated with either the signal, the 
initial-condition noise of the system, or even with itself from 
one measurement to the next. 

FIG. 7 shows a series of Monte Carlo simulations of SNR 
using the dispersive -operator version of the method with 
uncorrelated instrument noise added to each measurement. 
The spin polarization is again fixed at P=10 -5 . The value of 
the scaling parameter t, is chosen to be one half of the 
root-mean-square sum of the quantum noise and instrument 
noise, which is not necessarily optimized. Each curve rep- 
resents a fixed ratio of instrument noise o instr to quantum 
noise a Iz , where the quantum noise is a function of N. 
Curves for o ins Ja Iz equal to 0.01, 0.03, and 0.1 overlap with 
one another as indicated by curve 710 and produce the 
highest SNR at every N. Curve 720 represents the SNR for 
a insJ a iz = 0-3. Curves 730, 740, and 750, in order of decreas- 
ing SNR, represent values of o ins Ja Iz of 1, 3, and 10, 
respectively. In comparison, dashed lines show the SNR for 
the simulation of the corresponding first-order experiment 
using the same values for instrument noise as used in the 
second-order calculations. SNR decreases from one dashed 
line to the next for increasing instrument noise relative to 
quantum noise. 

As an example of the advantage in applying the present 
method to improve the SNR in a realistic experiment, 
quantum-noise-limited spectra have been simulated for the 
multiplication version of the method and the corresponding 
conventional first-order method assuming a cubic sample of 
a silicate crystal 2 microns on a side with natural abundance 
29 Si (spin Vi). The number of spins used in the calculation is 
N=6xl0 9 and room temperature polarization is assumed in 
a static magnetic field of 2 Tesla, which corresponds to 
P=1.4xl0 -6 . A 512-point time domain transient is simulated 
for each spectrum, with 1000 repetitions of the experiment. 

FIGS. 8A and 8B show approximated quantum-noise- 
limited spectra for (a) first-order pointwise NMR and (b) 
second-order pointwise NMR. The simulation indicates a 
significant improvement in signal-to-noise ratio of the 
second-order correlation NMR measurement over the con- 
ventional first-order measurement. It is noted that as the 
number of spins decreases, the signal in the first-order 
experiments disappears quickly into the noise, whereas the 
quantum-limited SNR in the second-order experiment using 
techniques of the present invention remains substantially 
unchanged. 

The performance of the present measuring technique can 
be enhanced by binning different second-order correlation 
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functions, i.e., the products S^), prior to further analysis, 
according to the value of the first and/or second measure- 
ment. Such binning may become increasingly useful as the 
accuracy of the measurements approaches the accuracy 
needed to assign the spin system to particular eigenvalues 5 
(or narrow subspaces of possible eigenvalues) of the mea- 
sured observable. In general, qualitatively different infor- 
mation can be obtained by calculating from the same data 
different averages <S 2 (t 1 )> using different linear combina- 
tions of the products of random deviates so binned. 10 

For example, by adding together only those measure- 
ments in which a particular eigenvalue of l z is measured both 
before and after the evolution period, an autocorrelation 
function of the projection operator onto this manifold can be 
obtained. This autocorrelation function, which is obtained, 15 
e.g., by an evolution period pulse sequence similar to that of 
FIG. 3, can include the frequencies of magnetic-dipole 
forbidden coherences. Hence, multiple -quantum spectra of 
uncoupled spins becomes possible without either the 
extended preparation or the mixing period normally needed 20 
in time -domain Fourier transform spectroscopy. Therefore, 
the present measuring technique does not require a spin- 
coupling network and can in general apply equally to 
coupled or uncoupled spins. Also, because the dephasing 
rates limit the length of preparation and mixing periods a 25 
procedure that replaces them by before and after observa- 
tions can make practical multiple-quantum spectroscopy of 
spins which are too weakly coupled to be put into the desired 
multiple-quantum superposition by way of their couplings. 
These weak couplings, providing in well-known ways infor- 30 
mation on molecular structure, would be apparent in the 
second-order interferogram as a function of t ± . 

Various embodiments of the present technique can be 
applied to mitigate various problems in a range of measure- 
ments that can be repeated on an ensemble. One important 35 
source of noise that is addressed is the uncertainty in 
measured quantities of interest that results from the distri- 
bution of initial conditions of the ensemble of systems under 
test. Many of the measurement problems addressed are most 
commonly encountered when the system under test or its 40 
component parts are microscopic in size. Of particular 
interest is spectroscopy, the measurement of the internal 
states of matter by way of their interaction with electromag- 
netic fields. 

For example, the present technique can be applied to various 45 
measuring schemes in coherent spectroscopy. Coherent 
spectroscopy, including time-domain coherent spectroscopy, 
provides a useful analyzing tool for many applications. See, 
Abragam, Principles of Nuclear Magnetism. Oxford, Clar- 
endon Press (1961) and Steinfeld, Molecules and Radiation: 50 
An Introduction to Modern Molecular Spectroscopy. 
Cambridge, Mass., MIT Press (1985). The term coherent 
applies to the nature of the time evolution of the system 
under test and/or reference systems within the spectrometer 
with which the quantities to be measured interfere. The 55 
property of the system, on which each measurement 
depends, may be a signed quantity or amplitude, such as an 
electromagnetic moment or mechanical displacement, or a 
positive quantity, such as an intensity or power or rate of 
photon counting. We note that in general these valuable 60 
methods are first order methods and could in a manner 
apparent to practitioners of these methods be converted into 
second order methods under the current invention. Such 
conversion would allow the goals of such methods to be 
obtained with sensitivity advantages of the present methods. 65 
The present technique may be used to incorporate any 
spectroscopic experiment by arranging that the desired evo- 
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lution of the spectroscopic system occurs during t 2 (i.e., 
between two measurements on the system). Preferably, the 
present technique is used in a time-domain experiment in 
which the interval t 2 is incremented and the spectrum is 
obtained by Fourier transformation with respect to t 1 . Mul- 
tidimensional experiments and frequency-domain experi- 
ments and any other time-domain experiments may also be 
incorporated in this period in ways that are obvious exten- 
sions. In some cases, other sets of experimental variables, 
such as time intervals or the frequencies, amplitudes, or 
phases of irradiation, will substitute for t 2 in ways that are 
known to practitioners of modern spectroscopy. In all these 
cases of the present method which can be described as 
converting prior art first-order measurements to second - 
order measurements, the algorithm is to make the before 
measurement in the first detection period, proceed with the 
steps of the first-order prior art (including its detection 
period, which comes the after measurement), and construct 
the possible second-order correlation functions of interest by 
repetition, as prescribed by the first-order procedure being so 
adapted. 

The present technique is generally applicable to any 
method of measuring spin magnetization or other internal 
properties of molecules or other systems. It is of value 
whenever the fluctuations of the observable of interest are 
comparable or larger than its expectation value, so pointwise 
first-order correlation function measurements or real-time 
methods would suffer from poor reproducibility due to the 
uncertainty of the initial condition. This technique does not 
require that the evolution of interest occurs during detection. 
This would, in general, prevent the separate optimization of 
sensitivity, resolution and information content; however, 
these optimizations are readily achieved within the context 
of the present method. 

The generalization of this procedure to involve more than 
two measurements separated by some manipulation of the 
spectroscopic system is also useful. For example, if in the 
magnetic resonance embodiment the “before” measurement 
resulted in an estimate [l z ] b whose magnitude was much 
smaller than the uncertainty in the total spin angular 
momentum, then a pulse (e.g., a jt/2 pulse) or other manipu- 
lation or evolution could be used to change the measurement 
statistics and another “before” measurement taken to replace 
the first “before” measurement, with the possibility that the 
magnitude of this second trial is larger and would thus 
increase the signal-to-noise ratio in the subsequent correla- 
tion function with the after measurement. Such a strategy, 
and many related ones in which other known types of 
spectroscopic manipulation would be interspersed with 
measurements, can be used for measuring correlation func- 
tions of order higher than two. 

Another context in which higher-order correlation func- 
tions are used is the analysis of the nonlinear response of 
matter to light to up to several incident pulses of light. Here 
the measurement of the emitted light takes place during or 
after the time evolution of interest. Elements of the method 
of the present invention which are absent from that work are 
repeated measurements interspersed with coherent manipu- 
lation of the system on the time scale of its relaxation and the 
multiplication of quantities derived from such separate mea- 
surements in constructing the spectra or time-domain inter- 
ferogram of interest. 

As illustrated in the above examples, an important area of 
application of the present invention will be magnetic reso- 
nance. Second and higher order correlation functions have 
been used in NMR measurements. In the method of sto- 
chastic magnetic resonance, a pseudo-random excitation of 
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the spin system near its resonant (Larmor) frequencies elicits 
a pseudo-random response in the transverse (precessing) 
magnetization of the spins. The response as measured and 
interpreted is proportional to the mean polarization of the 
ensemble at the start of the pulse sequence. The randomness 
is only apparent, and is in any case induced by the well- 
defined coherent incident irradiation. It is a determinate 
response of the system; it is not noise. 

One of the noise sources in NMR measurements is spin 
noise, the intrinsic quantum statistical uncertainty in the 
measurement of an observable of a spin system (e.g., its 
magnetization). The theoretical existence of spin noise in 
magnetic resonance measurements is well known. In its 
simplest guise it is a microscopic analog of the statistics of 
flipping coins, with the two orthogonal quantum states of a 
spin \=Vi playing the role of heads and tails. The mean 
polarization P of the ensemble (-1<P<1) plays the role of 
bias in the coin, with P=0 being the unbiased case. It is the 
mean value observed in first-order measurements of the 
magnetization normalized by dividing through by the spin 
quantum number I, proportional to the magnetic moment per 
spin, and the number of spins N. The quantities P and N fully 
specify the most accessible type of spin ensemble, that 
described as independent spin paramagnetism, which is a 
good approximation to the spin ordering obtained by equili- 
bration when the interaction of each spin with the a magnetic 
field is much stronger than the interactions between spins 
and at temperatures T where kT (with k the Boltzmann 
constant) is much greater than the energy of these interac- 
tions. The variance in P is the spin noise and has been 
recognized as a theoretical limit to the sensitivity of mag- 
netic resonance, although it is not clear that it has ever yet 
been the dominant noise in any practical attempt at obtaining 
a magnetic resonance spectrum. There have been some 
theoretical efforts to imagine ways of preparing an ensemble 
of spins in such a way that a first order measurement of P 
would have a smaller variance than that for an ensemble 
characterized by independent spin paramagnetism. 

Many known magnetic resonance spectroscopy experi- 
ments use first-order measurements of such spin polarization 
as a function of time or of the frequency of perturbing 
irradiation. The most notable exception is the experimental 
demonstration of spin noise in steady-state observations of 
rf magnetization. In these experiments, the zero-field spec- 
trum was obtained as a contribution to the noise in the 
SQUID-detected flux through a pickup coil near the sample. 
This spectrum was thus a power spectrum of the spontane- 
ous fluctuations of the magnetization with no applied radia- 
tion and thus is an example, in the frequency domain, of a 
second-order correlation function. 

A reported observation of the ESR of a single spin as a 
peak in the power spectrum of the radio-frequency current 
between a surface and an STM tip may be another analogous 
example. See, Manassen et al., Phys. Rev. Lett. 62 (1989) p. 
2531. The sensitivity for this steady -state, second-order 
approach has been quantitatively analyzed for detection by 
magnetic induction. Sleator et al. (1987). “Nuclear-Spin 
Noise and Spontaneous Emission.” Physical Review B Con- 
densed Matter 36: 1969-1980. The signal is the power 
spectrum of the fluctuations of that rf component of mag- 
netization which couples to the detector. For the hypotheti- 
cal case of a system with a sufficiently long spin-lattice 
relaxation time it has been argued that such an essentially 
noninvasive steady-state observation could be the preferred 
method of observation. It may not have been the preferred 
method in any realized situation. 

The present technique also uses the spin noise as part of 
the signal in the limit where P is negligible relative to its own 
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fluctuations. However, the present invention has no appli- 
cation to the experimental paradigm of measuring a spec- 
trum solely by way of the contribution of the steady-state 
spin noise to the frequency domain power spectrum of the 
5 observable. Rather it is aimed at the wider class of time- 
domain spectroscopic methods. It avoids steady-state by 
design and is indeed most useful when an entire repetition of 
the experiment can be performed in a time short compared 
to that at which the spin system reaches steady state by 
10 virtue of its interaction with the environment and detector. 

The present technique can be applied to measurements in 
force -detected magnetic resonance (FDMR). Conventional 
force -detected NMR techniques require an entanglement 
between the trajectories of individual molecules and their 
15 internal states throughout the experiment. The interaction 
between a laboratory field gradient and a permanent or 
transition dipole of the molecule of interest is used to encode 
spectroscopic information on the molecular trajectory, 
which is then measured. The ability of these conventional 
20 methods to provide spectroscopic resolution depends on the 
essentially collisionless nature of these trajectories. This 
usually limits their use to molecular beams and dilute 
trapped phases. In addition, the spectral resolution is often 
limited by time-of -flight through the apparatus or by the 
25 molecular collision time. This is in contrast to the more 
desirable situation of a high resolution-spectroscopy in 
which the achievable spectral linewidth is limited only by 
the lifetime or dephasing of the superpositions of internal 
states. In common with second-order methods, these force - 
30 detected methods, acting on single molecules with entangled 
trajectories, do not require any population differences (e.g., 
mean spin polarization P) between the states connected by 
the resonant irradiation. The entanglement allows this to be 
achieved despite there being a single measured quantity in 
35 these examples, for example the flux of molecules at a 
detector which terminates the trajectory. The present tech- 
nique does not require such entanglement and therefore can 
be advantageously used to measure NMR signals. Force - 
detected NMR methods can be used for ion spectroscopy in 
40 which the measurement of a trapped ion frequency before 
and after a period of spectroscopic evolution provides the 
information necessary to determine spectroscopic transition 
probabilities in an intermediate evolution period. See, 
Pizarro and Weitekamp (1992). Bull Magn. Reson. 14: 220 
45 and U.S. Pat. No. 4,982,088. These are second-order time- 
domain spectroscopic experiments which require that the 
members of the spectroscopic ensemble (e.g., molecules) be 
distinguished from one another during any part of the 
experiment by measurably different trajectories. The present 
50 technique can be advantageously applied in such ion spec- 
troscopic measurements for measuring states of matter 
where the molecules are rapidly colliding (e.g., fluids) or are 
strongly coupled in their motions (e.g., solids). This is 
because it is not required to have any degree of freedom with 
55 which the target observable need remain entangled over that 
period of the experiment which determines the spectroscopic 
resolution. 

The present technique may also be applied to forced- 
detected NMR measurements in which a sample is placed in 
60 a magnetic gradient field rather than a homogeneous field as 
in the previous examples. Such NMR systems and methods 
are disclosed by, for example, Sidles et al., “Magnetic- 
Resonance Force Microscopy.” Reviews Of Modern Physics 
67: 249-265(1995); Rugar et al., “Force Detection Of 
65 Nuclear-Magnetic-Resonance.” Science 264: 1560-1563 
(1994). A common feature is that the spin magnetization is 
modulated so as to drive a mechanical harmonic oscillator to 
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which it is coupled through magnetic forces. In the method 
of magnetic resonance force microscopy (MRFM). See, 
Sidles, 1995. The spin density of a sample in a large 
magnetic field gradient is imaged by successively bringing 
different regions (sensitive slices) of the sample into 
Larmor-frequency resonance with an rf magnetic field which 
is repeatedly swept through resonance so as to invert the 
longitudinal magnetization near the audio frequency of the 
harmonic oscillator. Either the sample or the ferromagnetic 
particle providing the field gradient can be the harmonically 
moving part. This method has been implemented by using 
conventional first-order correlation measurements, nomi- 
nally measuring the mean net polarization NP of each 
sensitive slice. It has been proposed that this method is 
extensible to the imaging of individual electron and nuclear 
spins, and it was recognized that such a measurement would 
not depend on the average spin polarization. Many of the 
instrumental requirements for such experiments have been 
detailed, and quantitative estimates of the anticipated signal- 
to-noise ratio (SNR) have been made that do not include spin 
noise and are dominated by Brownian motion of the 
mechanical oscillator. See, Sidles, J. A. and D. Rugar (1993). 
“Signal- to -Noise Ratios In Inductive and Mechanical Detec- 
tion Of Magnetic-Resonance.” Physical Review Letters 70: 
3506-3509. Omission of the spin noise contribution to the 
total noise limits the applicability of such SNR estimates to 
the case where exactly one spin (or none) interacts substan- 
tially with the measurement apparatus. 

Such limitations can be mitigated by using the present 
technique. Setting aside the difficulties of arranging for this 
to be the case, if it were achieved then the MRFM experi- 
ment would provide one bit of information on each repeti- 
tion; either the spin is present in the sensitive slice or it is 
not. This type of experiment can be viewed as a special case 
of a second-order correlation function in which the two 
factors are a single measurement whose value is squared. 
Were such prior art methods of MRFM to be applied to 
sensitive slices with P(N) 1/2 «1, but N >1, spin noise would 
dominate instrumental noise and previous (Brownian) SNR 
estimates would not apply. The present invention includes 
strategies for making more practical the generalization of 
MRFM to such sensitive slices containing arbitrary numbers 
of spins. Spin noise is an important consideration in the 
optimization of the method and in the resulting SNR. A 
qualitatively different analysis is needed and provided 
herein, indicating the limitations of the prior art MRFM 
methods and conception. The analysis also includes second- 
order adaptations of imaging methods such as Fourier zeug- 
matography and back-projection, which are advantageous, 
but incompatible with MRFM apparatus. 

A major limitation of MRFM apparatus as a magnetic 
resonance spectrometer is that there is a large field gradient 
at the sample. Magnetic resonance spectroscopy may be 
contrasted with magnetic imaging methods as being those 
experiments in which the frequencies of spin evolution to be 
determined are dominated by parameters (chemical shifts, 
spin-spin couplings, quadrupole interactions) other than the 
location of the spin along the field gradient. Indeed the vast 
majority of magnetic resonance methods require the mini- 
mization of the field gradient during all or part of the 
experiment. Although low resolution solid-state spectro- 
scopic information has been obtained with MRFM 
apparatus, this necessitates a magnetic field gradient far 
below that which is optimal for sensitivity. In such an 
apparatus, sensitivity is proportional to the static field gra- 
dient across the sample. A more extended discussion of the 
disadvantages of the static field gradient of the MRFM 
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method has been given. An additional problem worth noting 
is that the presence of the field gradient makes the resonance 
frequency such a strong function of molecular position that 
applications to fluids would suffer from severe loss of signal 
5 related to the short residence time of the spin-bearing 
molecule in the sensitive slice of a gradient large enough to 
offer a sensitivity advantage over inductive detection meth- 
ods. 

The BOOMERANG method of force-detected NMR uses 
10 a homogeneous static field at the sample to obviate the 
problems in other force-detected NMR methods and is 
especially advantageous in achieving higher SNR for small 
samples. In the BOOMERANG apparatus and method, the 
sensitivity and resolution are essentially independent, and 
15 both can generally be made better by orders of magnitude 
than in MRFM. A combination of the BOOMERANG 
method and the present technique can be used to fully utilize 
the advantages of both BOOMERANG and the present 
technique. 

20 FIGS. 9A and 9B are flow charts for an operation process 
in implementing the present technique using the BOOMER- 
ANG scheme for force -detection of NMR. The process may 
be divided into two parts: (1) preparation of sample and 
detection hardware; and (2) performing measurements and 
25 processing data using the present technique. 

FIG. 9 A depicts the operation steps in preparing a BOO- 
MERANG experiment. Certain details of these steps are 
disclosed in the incorporated references, U.S. patent appli- 
cation Ser. No. 08/872,528 and Leskowitz et al., “Force- 
30 detected magnetic resonance without field gradients,” Solid- 
State Nuclear Magnetic Resonance, Vol. 11(1,2), pp. 73-86 
(1998). For a given sample, it may be desirable to divide a 
sample into multiple pieces and use different detectors to 
measure different pieces. This can improve the detection 
35 sensitivity and increase the SNR. This is shown by step 910. 
Next, in steps 912-918, a composite magnet array, one or 
more oscillation sensors and other components are designed 
for a BOOMERANG apparatus to accommodate the sample. 

FIG. 9B is a flowchart for carrying out measurements and 
40 processing data using the present technique based on the 
BOOMERANG apparatus to obtain NMR spectra, NMR 
images, quantum computations, or other quantities of inter- 
est. At step 920, a “before” measurement is performed to 
measure the first response associated with the longitudinal 
45 magnetization of the sample. At step 921, the magnetization 
of the sample is determined. Next at step 923, an NMR pulse 
sequence is applied to the sample to allow the spins to 
evolve for obtaining spectroscopic measurements, NMR 
imaging, quantum computing and other effects. It may be 
50 necessary to further prepare the sample prior to the step 922. 
Such additional processing step is illustrated by the step 922 
in which one or more “before” measurements may be 
performed. After step 923, an “after” measurement is per- 
formed to obtain a second response from the sample on the 
55 longitudinal magnetization (step 924). If additional data 
points are needed, step 925 is carried out to allow the sample 
to recover in a “waiting period” and then steps 920-924 are 
repeated. Otherwise, at step 926, the “before” and “after” 
measurements are transformed to form desired correlation 
60 functions which are Fourier transformed for data extraction. 
More specifically, the measured time-domain transients are 
encoded pointwise by taking time-correlated measurements 
with a varying time t 1 in the evolution period (i.e., the third 
epoch in FIG. 1A) of the experiment and Fourier- 
65 transforming with respect to t 1 . 

Another experimental choice, which takes on special 
importance in the case where the initial condition noise is 
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dominant, is the fraction of possible observables which are 
included in a particular observation. In first-order spectro- 
scopic experiments, the SNR is ideally not affected by 
whether the observation made includes many spectral lines 
or only those of actual interest. This is because the noise 
relevant to a particular line is that which overlaps it in 
frequency and under typical conditions this is either additive 
in nature (e.g. thermal noise in the detector) or is associated 
with the transition (e.g. photon counting uncertainty intrinsic 
to that line). We define pointwise methods as those in which 
the spectrum is obtained as a Fourier transform with respect 
to an evolution period time variable L which is not also the 
detection time variable. In pointwise methods, the absence 
of correlation in noise between successive points in t 1 
guarantees that the noise in this conjugate frequency vari- 
able is “white,” or frequency independent. Thus, if the noise 
includes uncertainty due to fluctuations of the signal itself, 
as distinct from additive noise from the detection process, 
then there is a motivation to limit the signal to that which is 
advantageous [Weitekamp, 1983 #4] and to suppress that 
signal which is of no interest, but which would otherwise 
serve to contribute additional white noise. This noise is 
known as t ± noise and is ascribable in NMR to instrumental 
instability. A similar strategy is at play when difference 
experiments are used to isolate the effect of some experi- 
mental manipulation and when solvent lines are suppressed. 
Methods that can be used to minimize this problem by 
limiting the signals to those of interest are filters, frequency 
selective pulses, nk-quantum selective, differential relax- 
ation times, and inhomogeneous magnetic fields arranged to 
dephase unwanted contributions. The present technique can 
be combined with any of these techniques to improve the 
SNR of the measurements. The difference is that the present 
technique uniquely addresses the problem of fluctuations of 
signal that lead to noise which is intrinsic to the system, as 
opposed to being amenable to amelioration by improving the 
stability of the experimental manipulations. 

Yet another application of the present technique is quan- 
tum computing. A general computer may be regarded as a 
device that transforms input information into output infor- 
mation. Other operations performed by physical computing 
devices, such as rendering screen displays or printed pages, 
may be viewed as side effects of the “computing” or 
information-transforming process, and these side effects are 
irrelevant for the present discussion. In an ordinary, “clas- 
sical” digital computer, information takes the form of a 
string of bits, each of which may assume one of the values 
0 or 1. In contrast, a quantum computer processes quantum 
information in the form of qubits (quantum bits), each of 
which may be in a continuous linear combination, or 
superposition, of two states that may be defined as 0 and 1. 
The correspondence with the quantum mechanical two -level 
system is particularly useful here, and physical manifesta- 
tions of simple prior-art quantum computers and proposals 
for quantum computers are based on interacting two-level 
systems (e. g., spins or pairs of atomic electronic states). 
See, Monroe et al., Phys. Rev. Lett. 75 (1995) p. 4714 and 
Lloyd, Science 261, (1993) 1569. 

Interest in quantum computers surged when discovered a 
quantum computing algorithm for factoring integers that 
operates in a time that scales as a polynomial in L, the 
number of bits of input. Shor, in Proceedings of the 35 th 
Annual Symposium on the Foundations of Computer 
Science, (IEEE Computer Society, Los Alamitos, 1992), p. 
124. The factoring of large integers is an area of intense 
research, largely because of its connection to cryptography, 
but the best known algorithms that can be run on classical 
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computers require computing times proportional to exp(L 17 
3). The power of quantum computing algorithms lies in their 
use of input states that are superpositions of the ordinary 0 
and 1 states of classical computers. A suitably designed 
5 quantum computation on a superpositional input state may 
be viewed as a parallel computation on all possible input 
states. This parallelism is a key feature of quantum 
algorithms, and it allows polynomial-time computation 
where the best known classical computations require expo- 
10 nential time. There are at present few known quantum 
algorithms. See, Shor, supra.; Grover, Phys. Rev. Lett. 79 
(1997) p 325; and Brassard, Science 275, (1997) p 627. 

Quantum computers may be implemented by using indi- 
vidual two -level atoms as the qubits. A useful development 
15 in the field has been the introduction of ensemble quantum 
computing, in which “pseudo-pure” states of macroscopic 
systems serve as the qubits. See, e.g., Gershenfeld et al., 
Phys. Comp. 96, Proc. Of the Fourth Workshop on Physics 
and Computation, 1996, p. 134. These states are prepara- 
20 tions of a macroscopic system’s density operator that obey 
the same dynamics as states of a single microscopic system. 
Pseudo -pure states of, for example, nuclear spins in liquids 
have some advantages over single-atom states, namely that 
they are more robust to decoherence and are experimentally 
25 easier to prepare and manipulate. Simple quantum gates 
have been demonstrated in multiple -pulse NMR experi- 
ments with liquid samples at room temperature. Recent 
work has shown that a “universal quantum computer” is 
reducible in theory to suitably designed NMR experiments 
30 on suitably designed spin systems. Thus, one can in theory 
reduce doing arbitrary quantum computations to doing NMR 
experiments. 

Spin noise has been identified as a major challenge to be 
overcome in making ensemble quantum computing with 
35 NMR useful. See, Warren, “The Usefulness Of Nmr Quan- 
tum Computing.” Science 277: 1688-1689(1997). The main 
problem is that, in the regime where quantum computing 
might be superior to readily available classical computing, 
preparation of the required initial states is very inefficient. 
40 For example, in a computation requiring 100 qubits, the 
possible number of states for this system is 2 100 sl0 30 , and 
hence in a kilogram sample of ~10 24 molecules, a given 
initial state is usually unpopulated. This uncertainty in the 
initial condition is an aspect of the spin noise that is 
45 addressed by the present invention. By measuring the 
second-order correlation function of the magnetization, use 
is made of an arbitrary initial condition, and each shot of the 
experiment contributes to the signal. Many published 
ensemble quantum computing methods are described in 
50 terms of the first-order NMR signal. 

The present technique has been shown to enable the 
conversion of arbitrary multiple -pulse NMR experiments 
from first order experiments to second order experiments. 
Here we include all those NMR experiments for quantum 
55 computation. 

One procedure for converting quantum-computing NMR 
experiments to the second-order measurement methods of 
the present technique is as follows. NMR methods for 
quantum computing rely on detection of precessing trans- 
60 verse magnetization. The relevant quantum-mechanical 
expression for the first-order signal, in the representation in 
which the initial density operator (p(0)=l+bl^.) is taken to be 
that of the system after a preparatory jt/2 pulse, is 

65 S(f ; a)=<I + (f ; a)>=Tr{p(0)I + (f ; a)}=foTr{I JC I + (r ; a)} 

The parameter b represents the polarization of the spin 
system, and a represents all other parameters (usually time 
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periods and pulse phases) that characterize the evolution 
imposed on the system by the experimenter. To increase 
sensitivity, the same information that is available from this 
experiment may be obtained pointwise in a BOOMERANG 
experiment, where p(0)=l+bl^ and the observable is 1^: 

Sfl(^i,a)=<I z (r 1 ,a)>=Tr{p(0)I z (r 1 ,a)}=foTr{I z I z (r 1 ,a)}. 

Note that the “dimension” of the experiment is increased by 
one in this step. The experiment is now in a form suitable for 
measurement of the second-order correlation function. A 
measurement of l z is inserted prior to the period t 1? and the 
signal is defined as 

S 2 (f 1 ,a)=<I z (0)I z (f 1 ,a)>=Tr{p(0)I z I z (r 1 ,a)}=Tr{I z I z (r 1 ,a)}. 

Note that these expressions are identical to within a propor- 
tionality constant. The terms linear in b and of higher order 
in b have been dropped from the expression for the signal, 
for the method works even with unpolarized samples for 
which 

p(0)=l+M z ~l. 

These terms could be retained, as they were in the exact 
calculations illustrating the sensitivity of second-order 
BOOMERANG, but the principle value of the second-order 
measurement may be when the terms are small, which is the 
regime where the present invention can have a large sensi- 
tivity advantage over first-order methods. 

The measurement of the second-order correlation func- 
tion in this method allows measurement of \ z with approxi- 
mately unit SNR per root shot. For some quantum comput- 
ing applications, very high order multiple-quantum 
coherences may be the target observables, and there will 
necessarily be reduction in signal in the preparation of these 
coherences and in their conversion to observable magneti- 
zation. In any case, there is advantage to be gained by 
performing measurements simultaneously on as great a 
number of identically prepared small systems as is practical. 
Because of the small dependence of the second order SNR 
on the number of spins, the present method enables the 
complex spectroscopy needed for quantum computing to be 
done down to the single molecule level. Indeed, the limiting 
case is each target molecule being observed by multiple 
detectors. In principle, the best SNR would be obtained by 
coadding the signals from N molecules with the measure- 
ment periods no longer than needed to make instrument 
noise negligible with respect to spin noise. This ideal is 
approached by taking a macroscopic sample and dividing it 
among a large number of miniature spectrometers prepared 
perhaps by microfabrication techniques. The sensitivity will 
increase as the square root of the number of molecules. This 
procedure may also be used to advantage in non-quantum- 
computing NMR experiments where the sample may be so 
divided. 

While the above approach is a compelling strategy that is 
readily related to established first-order methods, it is not 
necessarily sufficient for a given spectroscopic or quantum 
computing application, particularly as the number of qubits 
becomes high. Here it is anticipated that binning the second- 
order products according to the values of the “before” and/or 
“after” measurements will also be useful. This allows prepa- 
ration and detection of high-order multiple -quantum coher- 
ence in a manner that is essentially independent of the 
couplings between spins. These superpositions of states play 
a key role in quantum computing algorithms, but traditional 
routes to their observation are impractical for the numbers of 
coupled spins that are needed for presently conceived quan- 
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turn computing algorithms. Note that the second-order con- 
cept is valuable when the system cannot for practical reasons 
be prepared in a pure state. When the number of distinguish- 
able spin states is as high as is envisioned in the quantum- 
5 computing literature, this will be the case even as P 
approaches unity. Thus it is expected that the present inven- 
tion will also be useful in situations of high, but imperfect 
spin polarization. 

A second measurement problem which can be substan- 
10 tially eliminated using techniques that are completely com- 
patible with the present method is that of the distortion of the 
desired information about the measured system by the “back 
action” of the measurement apparatus on that system. The 
principal problem caused by the interaction is a shortened 
15 decay time for a superposition of states whose Bohr fre- 
quency or intrinsic decoherence time is thereby obscured. A 
related problem is that of “frequency pulling”, where the 
observed frequencies are not that of the system alone, but of 
the system as coupled to the measurement apparatus. Such 
20 back action is well known both in the case of systems whose 
states span a finite Hilbert space (e.g., the phenomenon of 
“radiation damping” in NMR) and those which are well 
characterized by nominally infinite Hilbert spaces (e.g., 
harmonic oscillators). In a few cases, methods have been 
25 described which partially address this problem. In the case 
of NMR, this includes active suppression of the feedback of 
spin-induced rf fields in the sample inductor. In the field of 
trapped ion physics, measurement schemes have been dem- 
onstrated in which the interaction with the detector is 
30 subsequent to an evolution period during which the ion 
motion of interest is substantially undamped. See, Cornell et 
al., Phys. Rev. A 41 (1990) p. 312; and Phys. Rev. Lett. 63, 
(1989) 1674. A related issue is the considerable theoretical 
literature on “quantum non-demolition” (QND) measure- 
35 ments. Braginsky and Khalili, Quantum Measurement, 
Cambridge, Cambridge University Press (1992). This theory 
gives a recipe for finding that subset of observables whose 
repeated or continuous measurement can lead to accuracy 
which is not, in principle, limited by the Heisenberg uncer- 
40 tainty principle. 

The present method enables the adaptation of the above 
technique to second-order measurements and in particular to 
those which have the time line of FIG. 1A. In general, the 
present method allows for the absence of back action during 
45 the generalized evolution period. When such back action 
would otherwise be problematic, the method allows for the 
decoupling of the measurement apparatus from the system 
under test during this period. Such decoupling can be 
achieved by known means which can reduce back action to 
50 a negligible problem. These include altering the resonance 
frequency of the detector from that which would be opti- 
mum during measurement, broadening the resonance of the 
detector, or introducing shielding or physical distance or 
altered orientations between the sample and the detector. 

55 The present method also allows for the indirect second- 
order measurement of observables for which no instrument 
is practical or available. This includes the creation of 
so-called “forbidden” or disallowed spectroscopic 
superpositions, such as those involved in multiple-quantum 
60 coherence. The present method provides a second -order 
means for measuring their time evolution. In general, these 
superpositions will not be detectable in a second-order 
analysis of a continuously measured allowed observable. 
Rather, the time line of FIG. 1A is essential. It is simply 
65 necessary to provide, as part of the generalized evolution 
period, a preparation period at the beginning to convert the 
density operator following the “before” measurement to the 
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desired (forbidden) operator and a mixing period at the end 
(after a period of evolution) to convert the evolved operator 
to the observable relevant for the “after” measurement. 

We note that the prescription for QND observables 
includes (in principle, though perhaps not in known 
examples) forbidden QND observables. These can be 
observed by the same procedure, with a judicious choice of 
the QND observable (forbidden or not). 

Another application of the present invention is the char- 
acterization of systems of oscillators. Some physical 
examples such as trapped ions, phonons, the electromagnetic 
field, and microfabricated mechanical and electronic devices 
are normally described as harmonic oscillators. Varying 
degrees of anharmonicity are normally present and may be 
part of the information of interest. Examples with notable 
anharmonicity are vibrational states of molecules and elec- 
tronic states. 

Recall from the example of spin noise that initial condi- 
tion noise begins to dominate the ordinary first-order signal 
when the distribution of total spin angular momentum in the 
initial condition significantly overlaps with zero, leading to 
an increasing number of measurements in the signal average 
(ergodic ensemble average) whose sign is opposite to the 
sign of the average. The second-order multiplication scheme 
corrects for this sign fluctuation and converts repetitions of 
the experiment from destructive to constructive contributors 
to the signal average. 

An analogous initial condition for the quantum- 
mechanical harmonic oscillator is the thermal equilibrium 
density operator, whose average energy fluctuations are 
equal in magnitude (~kT) to the average energy, and whose 
phase distribution is totally random. In contrast to spin 
systems (for which the Hilbert space is bounded), with 
harmonic oscillators it is often practical to obtain a large 
first-order signal simply by coherently driving the oscillator 
from its thermal equilibrium initial condition to an ampli- 
tude much larger than the unamplified thermal fluctuations. 
The first-order signal is then measured and, perhaps, Fourier 
analyzed. However, in some applications it may be impos- 
sible or impractical to drive the initial condition to an 
amplitude significantly greater than the thermal distribution. 
The sample may be sufficiently remote from or inaccessible 
to a driving field, in which case the distribution of “before” 
measurements alone will have significant overlap with zero. 
In either case, the initial condition noise may dominate the 
signal, requiring the use of the proposed second-order tech- 
niques in order to achieve superior SNR. 

A similar situation is the case where one would like to 
characterize a system secretly. The minimally invasive 
second-order approach consists in this case of coupling to 
the system under test only during the measurement period 
and for only as long as needed for an acceptable measure- 
ment. This reduces the chance of the measurement being 
detected, while not compromising the spectral resolution of 
the measurement, which is determined by the length of the 
generalized evolution period, during which no measurement 
is made. Minimal perturbation of the target is desirable in 
espionage applications. 

Another related case is where exciting the system to large 
amplitudes would be undesirable or dangerous due to 
sample heating or where the system is easily driven into 
inelastic regions that would cause irreversible damage. Even 
where permanent damage would not be caused to the 
system, in some cases a large coherent excitation would 
drive the system into nonlinear response. A large- amplitude 
method would therefore degrade the measurement of “small- 
signal” properties of the system relative to a measurement 
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with very weak perturbations introduced in the coherent 
evolution period. This may be a prime concern in “nano- 
electronics” applications. In all of these cases, the optimal 
excitation, measurement, and/or evolution schemes may be 
5 incorporated into the present second-order method. 

A related situation, perhaps the most important applica- 
tion in the context of “minimally invasive measurements,” is 
for systems where the back-action from the measuring 
apparatus significantly changes the system of interest. In that 
10 case, it is desirable to provide for the system a period t 1 
during which the system evolves in the absence of interac- 
tion with the detection apparatus or with a thermal bath 
through the intermediary of the detection apparatus. The 
evolution thus proceeds in the system’s “natural state,” free 
15 of damping, “frequency-pulling,” or nonlinearities induced 
by the detector. The sensitivity of second-order measure- 
ments in which the detection apparatus is decoupled from 
the system during the evolution period may exceed that of 
analogous first-order procedures, and both the sensitivity 
20 and resolution may exceed that of continuous-measurement 
procedures analyzed with second-order correlation func- 
tions. Another example of the application of the present 
technique is the analysis of the thermal fluctuations of 
stationary systems. See, Nyquist, H. (1928). Physical 
25 Review 32: 110 and McCombie, C. W. (1953). Rep. Prog. 
Phys. 16: 266. By observing the system, ideally with neg- 
ligible perturbation by the measurement apparatus, and 
performing spectral analysis of various second-order corre- 
lation functions, it is possible to learn about the character- 
30 istic frequencies and relaxation (damping) times of the 
system. In most laboratory situations this is not the preferred 
approach because the fluctuations of a system in a stationary 
state are typically smaller than those that could be tran- 
siently induced by an intentional perturbation. While one 
35 could in principle determine that a piano is out of tune 
without touching it, by simply listening to its thermal 
fluctuations, it is preferable from the viewpoint of sensitivity 
to play it and spectrally analyze the resulting coherent 
response. This latter approach is an example of a first-order 
40 correlation function. First-order correlation functions pro- 
vide much of the same information about the system as the 
second-order correlation functions. They can be averaged 
over repetitions if desired. However, when coherent 
responses of such a system are comparable to or less than the 
45 fluctuations of the stationary state, measurement based on 
the second-order correlation functions of the present tech- 
nique becomes important. 

The present technique can be applied to obtain second and 
higher order correlation functions of the intensities of elec- 
50 tromagnetic fields. In a quantum mechanical description, the 
measured quantities are the positive -definite expectation 
values of the photon number operators of the quantized field, 
as distinct from the raising and lowering operators whose 
complex expectation values correspond to the field ampli- 
55 tudes. Operationally, the second order correlation function 
of a field has the form of <S 2 > with the quantities S fc (q fc , t b ) 
and S fl (q fl , t a ) being, for example, photocurrents in two 
detectors or the same detector at different times. Thus it is a 
measure of photon coincidences. A principal application of 
60 such observations is second-order interferometry, for 
example in astronomy. Another application of second and 
higher order correlation functions of photon number opera- 
tors is the characterization of nonclassical states of electro- 
magnetic fields. Another area of application is the optical 
65 spectroscopy of dilute molecules. Orrit et al., “Optical 
spectroscopy of single molecules in solids” in Progress in 
Optics , Vol. XXXV, pp. 61-144, (1996). 
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Although the present invention has been described in 
detail with reference to the preferred embodiments, various 
modifications and enhancements may be made without 
departing from the spirit and scope of the invention as 
defined by the following claims. 

What is claimed is: 

1. A method for suppressing measurement noise in a 
system having a memory which decays by a characteristic 
relaxation time, comprising: 

performing a first measurement on the system to produce 
a first response indicative of a first state of the system; 10 

performing a second measurement on the system after an 
evolution duration subsequent to said first measure- 
ment to produce a second response indicative of a 
second state of the system, wherein said evolution 
duration is selected so that the system retains a degree 15 
of memory of an aspect of said first measurement when 
said second measurement is performed; and 

forming a correlation function including a relation of a 
first function that has a relation with said first response 
and a second function that has a relation with said 20 
second response, 

wherein said correlation function is configured to repre- 
sent information indicative of the system and to reduce 
at least one component in said first and second 
responses that adversely affect extracting said informa- 25 
tion. 

2. A method as in claim 1, wherein said first function is 
identical to said first response data and said second function 
is identical to said second response data so that said corre- 
lation function is a product of said first response data and 3Q 
second response data. 

3. A method as in claim 1, wherein said first function is a 
transformation of said first response defined by 


MU) = 


Z h 

V — 

h z 


Zh 

z 2 +n 9 


35 


where is a scaling parameter and l z represents said first 
response. 40 

4. A method as in claim 1, further comprising applying 
perturbation to alter a state of the system during said 
evolution period. 


5. A method as in claim 4, wherein the system has 
particles with non-zero nuclear spins and said perturbation 
to the system includes two jt/2 electromagnetic pulses with 
opposite polarizations. 

6. A method as in claim 1, further comprising applying an 
initial perturbation to configure the system in a desired 
initial state. 

7. A method for performing a nuclear magnetic resonance 
measurement on a spin system having particles with non- 
zero nuclear spin moments and a characteristic relaxation 
time, sequentially comprising: 

preparing the spin system to be in an initial state; 

performing a first measurement on the spin system to 
produce a first response indicative of a first state of the 
spin system; 

allowing the spin system to evolve for a desired evolution 
period; 

performing a second measurement on the spin system to 
produce a second response indicative of a second state 
of the spin system, wherein a time interval between said 
first measurement and said second measurement is 
configured relative to the characteristic relaxation time 
in such a way that the spin system retains a degree of 
memory of an aspect of said initial state upon comple- 
tion of said second measurement; and 

forming a correlation function including a relationship of 
a first function that has a relation with said first 
response and a second function that has a relation with 
said second response to represent information indica- 
tive of a property of the spin system, 

wherein said correlation function is configured to reduce 
an effect of a fluctuation associated with said initial 
state. 

8. A method as in claim 7, wherein said preparing the spin 
system is configured to make said initial state a nonequilib- 
rium spin state. 

9. A method as in claim 7, wherein said preparing the spin 
system is configured to make the nuclear spins of the spin 
system unpolarized prior to said first measurement. 





